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(54) Method of determining the model of oxygen precipitation behaviour in a silicon 
monocrystalline wafer, method of determining a process for producing silicon 
monocrystalline wafers using said model, and recording medium carrying a program for 
determining the oxygen precipitation behaviour model in a silicon monocrystalline wafer 



(57) A method of determining oxygen precipitation 
behavior in a silicon monocrystal through use of a pro- 
grammed computer. According to this method, an initial 
oxygen concentration of a silicon monocrystal, an impu- 
rity concentration or resistivity of the silicon monocrys- 
tal, and conditions of heat treatment performed on the 
silicon monocrystal are input, and an amount of precip- 
itated oxygen and bulk defect density of the silicon 
monocrystal after the heat treatment are calculated 
based on the input data. The method enables quick, 
simple, and accurate determination of an amount of pre- 
cipitated oxygen and bulk defect density in silicon during 
or after heat treatment. 
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Description 

BACKGROUND OF THE INVENTION 



5 Field of the Invention: 



[0001 ] The present invention relates to a method of determining oxygen precipitation behavior in a silicon monocrys- 
tal, to a method of determining a process for producing silicon monocrystalline wafers, and to a recording medium on 
which is recorded a program for determining oxygen precipitation behavior in a silicon monocrystal. 

10 

Description of the Related Art: 

[0002] In general, silicon semiconductor wafers for use in fabrication of semiconductor devices such as ICs and LSIs 
are produced from a silicon monocrystal grown in accordance with a Czochralski method (CZ method). Since a silicon 
15 monocrystal contains supersaturated interstitial oxygen, the silicon wafers produced from the monocrystal also contain 
supersaturated interstitial oxygen. Therefore, if such a silicon wafer undergoes heat treatment in a process of fabricat- 
ing ICs or the like, interstitial oxygen precipitates in the form of silicon oxides, so that a large number of fine defects are 
generated inside the wafer. 

[0003] When such fine defects stemming from oxide precipitates are present within an internal region (bulk region) of 
20 a wafer, they desirably act as getter sites that capture heavy metal impurities and the like through so-called intrinsic get- 
tering (IG) . By contrast, when such fine defects are present within a semiconductor-device fabrication region in the 
vicinity of the surface of the wafer, the fine defects hinder operation of semiconductor devices, with the result that the 
characteristics of the devices degrade, and the production yield is affected directly and adversely. 
[0004] The above-described problem occurs regardless of whether semiconductor devices are fabricated on a mirror- 
25 polished wafer or an epitaxial wafer that is produced through deposition of monocrystalline silicon on a mirror-polished 
wafer. 

[0005] Therefore, when semiconductor devices are fabricated on silicon wafers produced in accordance with the CZ 
method, determining oxygen precipitation behavior within a silicon monocrystal as well as a process for producing sili- 
con monocrystalline wafers in accordance with the determined oxygen precipitation behavior are considerably impor- 
30 tant for obtaining, at a desired yield, semiconductor devices having desired characteristics. 

[0006] Thus, the intrinsic gettering (IG) capability of a silicon wafer must be known. In order to express such IG capa- 
bility, there have been used an amount of precipitated oxygen, which is a decrease in the concentration of interstitial 
oxygen occuning due to precipitation during a heat treatment process, or a density of bulk defect formed during the heat 
treatment process. The amount of precipitated oxygen and the bulk defect density are known to strongly depend on 
35 concentration of interstitial oxygen within a silicon wafer (initial oxygen concentration within a silicon monocrystal pro- 
duced in accordance with the CZ method) and conditions of heat treatment that the silicon wafer has undergone (includ- 
ing ail thermal histories of silicon, such as a thermal history during growth of the silicon monocrystal in accordance with 
the CZ method and that of heat treatment that the silicon wafer has undergone in a device-fabricating process). 
[00071 Accordingly, a proper initial value of oxygen concentration of silicon wafers to be used and conditions of heat 
40 treatment for the wafers must be actually investigated in order to ensure that the wafers will have a desired gettering 
capability during or after a specific device-fabricating process and will have a proper amount of precipitated oxygen and 
a bulk defect density to thereby increase device yield. Conventionally, a proper initial value of oxygen concentration of 
silicon wafers and conditions of heat treatment for the wafers have been determined as follows: A large number of 
wafers having different initial oxygen concentrations and condition of heat treatment performed are prepared and sub- 
45 jected to a heat treatment actually performed in a device-fabricating process to be used or to a heat treatment that sim- 
ulates the actual heat treatment, under various heat treatment conditions. Subsequently, the amount of precipitated 
oxygen and the bulk defect density of the wafers are measured in order to determine a proper initial value of oxygen 
concentration of silicon and heat treatment conditions. The conditions of production of a silicon monocrystal according 
to the CZ method and the conditions of heat treatment for silicon wafers are determined in order to achieve the deter- 
so mined initial oxygen concentration and the determined heat treatment conditions. 

[0008] However, when such a conventional method is employed, there must be prepared a large number of wafers 
having different initial oxygen concentrations and thermal histories, and heat treatment must be performed for a pro- 
longed period of time. Further, since heat treatment conditions vary depending on the type and maker of semiconductor 
devices to be fabricated, the amount of precipitated oxygen and bulk defect density of wafers must be examined through 
55 actual heat treatment whenever the type or maker of semiconductor devices changes. Therefore, considerable time and 
money are needed in order to determine an initial oxygen concentration and heat treatment conditions for obtaining a 
proper amount of precipitated oxygen and bulk defect density. 

[0009] In addition, since an initial oxygen concentration and heat treatment conditions have been determined on the 
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spur of the moment and empirically, accuracy of determination has been insufficient, so that in many cases walers sub- 
jected to an actual production process do not have desired characteristics, due to various factors (e.g., difference in 
crystal production apparatus and heat treatment apparatus) Therefore, in some cases, it has been difficult to reliably 
obtain a desired initial oxygen concentration and heat treatment conditions. 

5 

SUMMARY OF THE INVENTION 

[0010] The present invention has been accomplished to solve the above-mentioned problems, and an object of the 
invention is to provide a method which can simply and accurately determine an amount of precipitated oxygen and bulk 
w defect density of silicon during or after heat treatment, within a considerably short time, through numerical calculation 
which is performed through use a programmed computer and on the basis of an initial oxygen concentration, heat treat- 
ment conditions, and other factors. 

[001 1] Another object of the present invention is to provide a method of determining a process for producing silicon 
monocrystalline wafers, in which the dependency of the amount of precipitated oxygen and bulk defect density on the 

15 various factor is determined in accordance with the above-described method, and in which various factors such as an 
initial oxygen concentration within a silicon monocrystal and heat treatment conditions are determined such that wafers 
will have a desired amount of precipitated oxygen and bulk defect density during or after a specific process. 
[001 2] The inventors of the present invention found that oxygen precipitation behavior; i.e. , an amount of precipitated 
oxygen and bulk defect density, greatly depend on the initial oxygen concentration of a silicon monocrystal, conditions 

20 of heat treatment performed on the monocrystal, and the impurity concentration or resistivity of the monocrystal; and 
made it possible to simply and precisely determine oxygen precipitation behavior through numerical calculation by stor- 
ing the above-described dependency relationship in a computer in a form of a program and data. 
[0013] Specifically, to achieve the above objects, the present invention provides a method of determining oxygen pre- 
cipitation behavior in a silicon monocrystal through use of a programmed computer. According to this method, an initial 

25 oxygen concentration of a silicon monocrystal, an impurity concentration or resistivity of the silicon monocrystal, and 
conditions of heat treatment performed on the silicon monocrystal are input, and an amount of precipitated oxygen and 
bulk defect density of the silicon monocrystal after the heat treatment are calculated based on the input data. 
[0014] The above-described numerical calculation through use of a computer enables quick and precise determina- 
tion of oxygen precipitation behavior without performance of an experiment in which heat treatment and the like are per- 

30 formed on silicon wafers. Especially, in the present invention, since an impurity concentration or resistivity of a silicon 
monocrystal is taken into consideration, along with an initial oxygen concentration and heat treatment conditions, the 
accuracy of the numerical calculation is improved greatly. 

[001 5] Preferably, the computer is programmed based on a model in which oxygen precipitation is expressed as a 
function of the oxygen concentration of the silicon monocrystal and the temperature and time of the heat treatment, and 

35 free energy produced due to formation of oxide precipitates depends on the impurity concentration of the silicon monoc- 
rystal. In this case, since only a small difference is produced between values obtained through actual heat treatment 
and calculated values, the oxygen precipitation behavior within a silicon monocrystal can be determined accurately. 
[0016] More preferably, the oxygen precipitation behavior within a silicon monocrystal is determined under the 
assumption that the free energy depending on the impurity concentration is strain energy. 

40 [001 7] The present invention also provide a method of determining a process for producing silicon monocrystalline 
wafers. According to this method, the initial -oxygen-concentration dependency, impurity-concentration or resistivity 
dependency, and crystal -thermal-history dependency of the amount of precipitated oxygen and bulk defect density are 
determined in accordance with the above-described method of the present invention. Subsequently, an initial oxygen 
concentration, impurity concentration, resistivity, and thermal history of the silicon monocrystal are determined such 

45 that wafers will have a desired amount of precipitated oxygen and bulk defect density during or after a specific wafer 
process. 

[0018] As described above, once various oxygen concentration, impury concentoration or resisitivity, conditions of 
heat treatment and the like are input, the initial -oxygen-concentration dependency, impurity-concentration or resistivity 
dependency, and crystal-thermal-history dependency of the amount of precipitated oxygen and bulk defect density are 

so determined in accordance with the above-described method of the present invention, an initial oxygen concentration, 
impurity concentration, resistivity, and thermal history of the silicon monocrystal can be simply determined such that 
wafers will have a desired amount of precipitated oxygen and bulk defect density during or after a specific wafer proc- 
ess. Preferable results are obtained when a silicon monocrystal having a thus-determined initial oxygen concentration 
is produced in accordance with the CZ method and conditions of heat treatment performed in a monocrystal growth 

55 process or a subsequent process are determined in accordance with the calculated value. 

[0019] The present invention further provides a recording medium on which is stored a program of determining oxy- 
gen precipitation behavior in a silicon monocrystal through use of a computer. The program instructs the computer to 
input an initial oxygen concentration of a silicon monocrystal, an impurity concentration or resistivity of the silicon 
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monocrystal, conditions of heat treatment performed on the silicon monocrystal, and to calculate, based on the input 
data, an amount of precipitated oxygen and bulk defect density of the silicon monocrystal after the heat treatment 
[0020] Since the program prepared on the basis of a model for carrying out determining procedures according to the 
present invention can quickly and precisely determine oxygen precipitation behavior, rf such a program is recorded on 
5 a recording medium, the program can be input into a computer for use at any time and any place, thereby providing con- 
siderable convenience. 

[0021] According to the present invention, if an initial oxygen concentration, impurity concentration or resistivity, and 
heat treatment conditions are known, an amount of precipitated oxygen and a hulk defect density of silicon during or 
after heat treatment can be simply and precisely determined, within a considerably short time, through numerical cal- 

10 culation performed by means of a programmed computer. 

[0022] Further, once the dependency of the amount of precipitated oxygen and bulk defect density of a silicon monoc- 
rystal on various factors is determined, an initial oxygen concentration, impurity concentration or resistivity, and heat 
treatment conditions of the silicon monocrystal can be determined such that wafers will have a desired amount of pre- 
cipitated oxygen and bulk defect density during or after a specific process. 

is [0023] Accordingly, the present invention obviates necessity of actually performing heat treatment for a prolonged 
period of time on a large number of wafers having different initial oxygen concentrations and thermal histories. There- 
fore, there can be quickly determined, at low cost, an initial oxygen concentration, impurity concentration or resistivity, 
and heat treatment conditions for yielding silicon wafers having a proper amount of precipitated oxygen and bulk defect 
density. 

20 [0024] Since the determination is performed neither on the spur of the moment nor empirically and is therefore accu- 
rate, when wafers are actually processed in order to fabricate semiconductor devices, the semiconductor devices can 
be fabricated at an improved yield. 



BRIEF DESCRIPTION OF THE DRAWINGS 

25 

[0025] 

FIG. 1 is a diagram showing temperature variation in a heat treatment that was modeled on a device fabricating 
process and was used in an Example of the present invention; 
30 FIG. 2 is a graph showing the results of the Example (relationship between initial oxygen concentration and amount 
of precipitated oxygen); 

FIG. 3 is a graph showing the results of the Example (relationship between initial oxygen concentration and bulk 
defect density); 

FIG. 4 is a graph showing the results of a comparison between amount of precipitated oxygen measured through 
35 actual heat treatment and calculated amount of precipitated oxygen shown in FIG. 2; and 

FIG. 5 is a graph showing the results of a comparison between bulk defect densities measured through actual heat 
treatment and calculated bulk defect densities shown in FIG. 3. 



DESCRIPTION OF THE INVENTION AND PREFERRED EMBODIMENT 

40 

[0026] The present invention will next be described in detail. 

[0027] The inventors of the present invention carried out various studies in order to find a way for determining, through 
numerical calculation, a proper initial oxygen concentration and heat treatment conditions more easily than in the case 
of a conventional manner in which wafers of different initial oxygen concentrations, ranging from a low oxygen concen- 
45 tration to a high oxygen concentration, are actually heat-treated for a long time under various heat treatment conditions 
in order to empirically find a proper initial oxygen concentration and heat treatment conditions. The present invention 
was accomplished based on the results of such studies. 

[0028] The present invention enables quick, low cost, and accurate determination of an amount of precipitated oxygen 
and bulk defect density through numerical calculation, which is performed upon input of only an initial oxygen concen- 

so tration, impurity concentration or resistivity, and heat treatment conditions. 

[0029] Especially, the present invention is characterized in that the accuracy of determination has been greatly 
improved through employment of a numerical calculation program that uses a combination of a conventional model for 
oxygen precipitation in a silicon monocrystal and an additional model representing that free energy (more specifically, 
strain energy) generated due to formation of silicon oxides depends on an impurity concentration or resistivity of the sil- 

55 icon monocrystal. 

[0030] By virtue of the above-described method, the initial-oxygen-concentration dependency, impurity-concentration 
or resistivity dependency, and crystal-thermal-history dependency of the amount of precipitated oxygen and bulk defect 
density are determined in advance through substitution of various values. Thus, by use of the thus-determined various 
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dependencies, an initial oxygen concentration, impurity concentration, resistivity, and thermal history of the silicon 
monocrystal can be determined such that wafers will have a desired amount of precipitated oxygen and bulk defect den- 
sity during or after a specific wafer process. Therefore, silicon wafers having a desired gettering capability can be easily 
and reliably obtained if conditions for a process of producing a silicon monocrystal according to the CZ method, and 

5 conditions for subsequent wafer processing and device fabricating processes are determined such that the determined 
initial oxygen concentration, impurity concentration or resistivity, and heat treatment conditions are attained. 
[0031] The expression "impurity concentration or resistivity" is used because an impurity that causes a problem in 
relation to highly pure silicon wafers generally used in fabrication of ICs and the like is a dopant, and the concentration 
of such an impurity directly affects the resistivity. 

w [0032] An embodiment of the present invention will next be described with reference to details of a model based on 
which a program is prepared. However, the present invention is not limited thereto. 

[0033] Conventionally, oxygen precipitation behavior has been expressed through use of a Fokker-Planck expression 
reported by Schrems. Here, a description will be given of an exemplary model based on this expression. 
[0034] In this model, oxygen precipitation behavior is statistically described by a function f(r.t) representing size dis- 
15 tribution o1 precipitates, wherein r is a radius of an oxide precipitate under the assumption that the oxide precipitate is 
spherical, and t is a heat treatment time. The partial differentiation of f(r,t) with respect to time is described by the fol- 
lowing Expressions (1) and (2), where are collectively called the Fokker-Planck expression: 

= -§- r J(r.t) (1) 



25 



J = -B(rj)f r +A(r,t)f (2) 



where A (r, t) and B(r, t) satisfy the following Expression (3): 



30 dr 



where k is Boltzmann's constant, and T is an absolute temperature. 

The Gibbs's free energy AG(r, t) of a precipitate having a radius r is given as a sum of a volumetric energy and an inter- 
35 facial energy as shown in Expression (4) below: 



3v o C% q 



(4) 



40 



where C eq is a thermal equilibrium concentration of oxygen in silicon, a is an interfacial energy at the interface between 
a matrix and an oxide precipitate, and V 0 is half the volume of an Si0 2 molecule. In Expression (4) representing the free 
energy AG, for the sake of simplicity, effects of stresses of precipitates and intrinsic point defects are not considered. 
45 [0035] The density C 0 (t) of residual interstitial oxygen is calculated in accordance with the following Expression (5) 
based on the assumption that the overall density C T of oxygen atoms is maintained. 



50 



C o (0= C r - J (£L r 3 )f{r.t)dr 



(5) 



where r mjn is the radius of one Si0 2 molecule. 

[0036] A(r. t) is approximately equal to a growth rate dr/dt of precipitate and is expressed by Expression (6) below: 

55 

A(r.f)-g =^plC 0 (0-Ci'(O] (6) 
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where D 0 is a diffusion coefficient of oxygen within silicon. 

[0037] Based on the assumption aAG/dr = 0 , which is a thermodynamical equilibrium condition at the interface 
between the matrix and a precipitate, the oxygen concentration C 0 ,f (r) at the interface is expressed by Expression (7) 
below: 

C ;,V)=Cexp(?g- 0 ) (7) 

[0038] In order to describe an initial condition and boundary conditions of the Fokker-Planck expression, a quasiequi- 
librium distribution function f 0 (r, t) is introduced. Through use of this function, the initial condition can be described as 
f(r, 0) » f 0 (r, 0)8rr mjn where 6 is a delta function. The boundary conditions include f(r min , t) = f 0 (r mjn , t) and 
f(oo, t) = 0. f 0 (r, t) is given by Expression (8) below: 

Wf.0-^C ( (0P^H^) (8) 



[0039] The above is the model based on the Fokker-Planck expression reported by Schrems (see M. Schrems, in oxv- 
20 oen in Silicon. Semiconductors and Semimetals. Chap. 10, Vol. 42, R. K. Willardson, A. C. Beer, and E. R. Weber, 
Series Editors, F. Shimura, Volume Editor, p. 391 , Academic Press, New York (1994)). 

[0040] However, in the Schrems's model, the coefficient P in Expression (8) is considered a variable that is affected 
by a thermal history during a crystal pulling process and by impurities other than oxygen, and therefore is handled as a 
fitting parameter used in order to match a calculation result to an experimental result. Therefore, an experiment for 
25 obtaining the coefficient P must be performed, after all. In addition, since a value of the coefficient P obtained through 
such experiment is valid only for certain conditions, the coefficient P must be obtained whenever the conditions are 
changed, which is still cumbersome. Moreover, the obtained results tend to involve a large error and therefore be inac- 
curate. 

[0041 ] In view of the foregoing, the inventors of the present invention carried out various studies and found that the 
30 results of various experiments can be expressed without adjustment by a fitting parameter, when the coefficient P is 
expressed as a function of a heat treatment temperature T and an interstitial oxygen concentration C 0 (t), as expressed 
by Expressions (9) and (10) below: 



35 



where a = 1.92 eV, 



P = P 0 exp(^) (9) 



Pomb{ 2M.xf no) 



where b = 2.97x10- 18 and c=5.45x10 17 cm* 3 when C 0 (t)>7x10 17 cm" 3 , and P 0 = 1.46x10 -19 when C 0 (t)^7x10 17 cm -3 . 
45 (see H. Takeno, T. Otogawa and Y. Kitagawara, J. Electrochem. Soc. 144, 4340 (1997)) 

[0042] As described above, since oxygen precipitation is expressed as a function of an oxygen concentration, a heat 

treatment temperature, and a heat treatment time, oxygen precipitation behavior can be determined without use of a 

parameter, such as the above-described fitting parameter, that must be obtained through an experiment. 

[0043] Although in some cases values obtained through use of the above-described method agree well with values 
so obtained through actual heat treatment, the values obtained through use of the above-described method sometimes 

involve considerable errors. Therefore, the above-described method is insufficiently accurate in determining oxygen 

precipitation behavior that precisely agrees with the actual behavior. 

[0044] In view of the foregoing, the inventors of the present invention investigated a cause of such errors and found 
that the concentration of impurities (dopant) in a silicon monocrystal greatly affects the oxygen precipitation behavior. 
55 That is, such errors are introduced when the above-described method is applied to wafers, such as wafers for epitaxial 
growth, in which impurities are doped at a high concentration. On the basis of this finding, the above-described model 
is modified in order to take into consideration an impurity-concentration dependency (i.e., in order to obtain a model in 
which an impurity-concentration dependency of free energy is taken into consideration). The inventors of the present 
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invention succeeded in creating a program for determining oxygen precipitation behavior that agrees well with that 

obtained through actual heat treatment. Next, a description will be given of the thus-created program. 

[0045] Expression (1 1) below is widely accepted in general for oxygen precipitation reaction in a silicon crystal: 

s (1 + a)Si = 20, o pl/= SiO s + at + Stress (11) 

where Si is a silicon atom, Oj is an interstitial oxygen atom, V is a void, Si0 2 is a silicon oxide molecule, and I is an inter- 
stitial silicon atom, a is the ratio of discharged interstitial Si atoms for every Si0 2 molecule, and 3 is the ratio of absorbed 

10 voids for every Si0 2 molecule. 

[0046] Since the volume of an Si0 2 molecule is about 2.2 times that of an Si atom, stresses generated due to forma- 
tion of Si0 2 are relaxed through absorption of voids and discharge of interstitial Si atoms. In a boron-doped crystal, 
since the volume of a boron atom is smaller than that of an Si atom, boron atoms are conceivably absorbed in Si0 2 in 
order to relax the stresses. From this idea, the inventors of the present invention derived the following Expression (12), 

15 in place of Expression (11 ), for a boron-doped crystal. 

(1 + a)S/=20,- + p\/+ yB = SiO 2 + at + Stress (12) 

20 where B is a boron atom, and y is the ratio of absorbed boron atoms for every Si0 2 molecule. 

[0047] With respect to the reaction of Expression (12), the linear misfit 0 of Si0 2 is expressed by Expression (13) 
below: 



25 



35 



e={ [i + (a + p)]£ 2 Sl+ rn e )3 - 1 (13) 

where 

£lsio 2 

is the volume of an Si0 2 molecule, Qsi is the volume of an Si atom, and Qb is the volume of a boron atom. 
[0048] Since 



is about 2.4 times fiB, the linear misfit 6 decreases with absorption of boron atoms, as well as with absorption of voids 
40 and discharge of interstitial Si atoms. 

[0049] Here, the stress energy AG 3 (strain energy) of an oxide precipitate having a radius r is expressed by Expres- 
sion (14) below: 

where p. is the shear coefficient of silicon, and K is the compression ratio of Si0 2 . 

[0050] According to Expressions (13) and (14), absorption of boron atoms into oxide precipitates causes a reduction 
so in the stress energy of the precipitates. 

[0051 ] Also, the ratio y of absorbed boron atoms is assumed to be proportional to the density [B] of boron doped into 
a silicon crystal as shown by Expression (15) below: 



55 



Y = d[B] (15) 

where d is a constant. 

[0052] In order to calculate oxygen precipitation behavior in a boron-doped crystal, the above-described model is 
combined with the model based on the Fokker-Planck expression. That is, when the strain energy of Expression (14) is 
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added to the free energy of Expression (4). the strain energy that depends on the impurity concentration is corrected, 
so that calculated values agree well with experimental values even in wafers in which impurities are doped at a high 
concentration. 

[0053] In an actual program, based on Expression (3) and Expression (4), to the right side of which is added Expres- 
sion (14), calculation is performed for Expressions (1) and (2) in accordance with a differential calculus in order to obtain 
the size distribution function f(r, t) Further, an amount of precipitated oxygen and bulk defect density are calculated from 
the thus-obtained function (r, t). 



10 



EXAMPLES 

[0054] The present invention will next be described by way of example, which should not be construed as limiting the 
invention. 



Example: 

15 

[0055] First, boron-doped silicon wafers having respective resistivities of 10, 0.02, and 0.01 n- cm were sliced from 
crystal ingots prepared in accordance with the CZ method and were then subjected to epitaxial growth. Subsequently, 
the relationship among initial oxygen concentration, amount of precipitated oxygen, and bulk defect density was found 
or determined through calculation by use of a programmed computer of the present invention, for the case where there 
20 was performed a heat treatment as shown in FIG. 1 that was modeled on a device fabricating process. 

[0056] FIGS. 2 and 3 show the results of the calculation. FIG. 2 shows the relationship between initial oxygen con- 
centration and amount of precipitated oxygen, while FIG. 3 shows the relationship between initial oxygen concentration 
and bulk defect density. 

[0057] In order to verify the calculation results of FIGS. 2 and 3, the silicon wafers were actually heat-treated, and 

25 their amount of precipitated oxygen and bulk defect densities were measured. 

[0058] Specifically, boron-doped silicon wafers having respective resistivities of 10, 0.02, and 0.01 n • cm were first 
sliced from crystal ingots prepared in accordance with the CZ method. The initial oxygen concentrations of the thus- 
prepared silicon wafers were measured in accordance with an infrared absorption spectroscopy (FT-IR) (for the wafer 
having a resistivity of 10 n • cm) or a gas fusion method (for the wafers having respective resistivities of 0.02 and 0.01 

30 n*cm). 

[0059] Subsequently, the wafers whose initial oxygen concentrations had been measured were subjected to epitaxial 
growth such that silicon monocrystalline film of about 5 urn thick was grown at 1 130°C. Subsequently, the heat treat- 
ment that had been modeled on the device fabricating process as shown in FIG. 1 was actually performed. 
[0060] After the heat treatment, the oxygen concentration (concentration of residual interstitial oxygen) of the wafer 

35 having a resistivity of 10 n« cm was measured in accordance with infrared absorption spectroscopy (FT-IR). Subse- 
quently, a difference between the measured oxygen concentration and the initial oxygen concentration was calculated 
in order to obtain an amount of precipitated oxygen. Further, the oxygen concentrations of the wafers having respective 
resistivities of 0.02 and 0.01 n • cm were measured in accordance with a method utilizing x-ray diffraction (see Japa- 
nese Patent Application Laid-Open (kokai) No. 8-75680). Further, the bulk defect densities of the wafers were meas- 

40 u red in accordance with an infrared laser scattering tomograph method. 

[0061] In the gas fusion method, silicon in bulk form is placed in a carbon crucible and heated in an inert gas atmos- 
phere to a temperature above the melting point of silicon, and the concentration of carbon monoxide or carbon dioxide 
produced from a reaction between oxygen within silicon and carbon of the crucible is detected by use of an infrared 
detector to thereby measure oxygen concentration. 

45 [0062] In the infrared laser scattering tomograph method, an infrared laser is radiated into a crystal and scattered light 
stemming from defects in the crystal is detected in order to detect the defects. 

[0063] The density of interstitial oxygen was measured in accordance with the standard of Japanese Electronic Indus- 
try Development Association (J El DA). 

[0064] FIGS. 4 and 5 show the results of the measurement. In FIG. 4 are plotted amount of precipitated oxygen meas- 
so ured through actual heat treatment and amount of precipitated oxygen obtained through calculation shown in FIG. 2. In 
FIG. 5 are plotted bulk defect densities measured through actual heat treatment and bulk defect densities obtained 
through calculation shown in FIG. 3. 

[0065] As is apparent from FIGS. 4 and 5, at all levels of resistivity, the amount of precipitated oxygen and bulk defect 
densities calculated in accordance with the present invention agree well with those obtained through measurement per- 
55 formed after actual heat treatment and the like. 

[0066] According to the present invention, if initial oxygen concentration, impurity concentration or resistivity, and heat 
treatment conditions are known, oxygen precipitation behavior can be determined. Therefore, once various oxygen con- 
centration, impurity concentration or resistivity, conditions of heat treatment and the like are input, the initial-oxygen- 



8 



EP 0 917 192 A1 



concentration dependency, impurity-concentration or resistivity dependency, and crystal-thermal-history dependency of 
the amount of precipitated oxygen and bulk defect density are determined in advance, an initial oxygen concentration, 
impurity concentration or resistivity, and heat treatment conditions of the silicon monocrystal can be determined such 
that wafers will have a desired amount of precipitated oxygen and bulk defect density during or after a specific wafer 
process. Therefore, silicon wafers having a desired gettering capability and a high production yield can be obtained if a 
silicon monocrystal having the thus-determined initial oxygen concentration and an impurity concentration or resistivity 
is produced according to the CZ method, and the process for producing silicon monocrystalline wafers is determined 
such that conditions for the growth of the monocrystal and subsequent heat treatment agree with the calculation results. 
[0067] If a program for causing a computer to determine oxygen precipitation behavior in a silicon monocrystal is 
recorded on a recording medium, the program can be input into the computer for use at any time and any place, thereby 
providing considerable convenience. For example, such a program is input to a computer that controls a crystal produc- 
tion apparatus in accordance with the CZ method or a computer that controls various types of heat treatment furnaces, 
thereby making the program considerably useful. 

[0068] The present invention is not limited to the above-described embodiment. The above-described embodiment is 
a mere example, and those having the substantially same structure as that described in the appended claims and pro- 
viding the similar action and effects are included in the scope of the present invention. 

[0069] For example, although the model of oxygen precipitation behavior used in the above-described embodiment 
is based on the Fokker-Planck expression, the present invention is not limited thereto; no limitation is imposed on the 
basic model insofar as the model allows determination of oxygen precipitation behavior from an initial oxygen concen- 
tration, impurity concentration or resistivity, and heat treatment conditions of silicon wafers. 

[0070] Further, although in the above-described embodiment the free energy that depends on impurity concentration 
is assumed to be strain energy, the present invention is not limited thereto; chemical potential and/or interfacial energy 
may also be taken into consideration. 

[0071] The expression "heat treatment conditions" used in relation to the present invention encompasses conditions 
of all heat treatments performed in a process for producing a silicon monocrystal according to the CZ method, a process 
for producing wafers from the crystal, and a process for fabricating semiconductor devices on the wafers. Therefore, the 
expression "determination of oxygen precipitation behavior in silicon wafers and a silicon monocrystal" is not limited to 
determination of oxygen precipitation behavior after a specif ic device fabricating process, and also encompasses deter- 
mination of oxygen precipitation behavior after crystal growth according to the CZ method, after wafer processing, and 
during the device fabricating process. 

Claims 

1 . A method of determining oxygen precipitation behavior in a silicon monocrystal through use of a programmed com- 
puter, characterized by comprising the steps of: 

inputting an initial oxygen concentration of a silicon monocrystal, an impurity concentration or resistivity of the 
silicon monocrystal, and conditions of heat treatment performed on the silicon monocrystal; and 
calculating an amount of precipitated oxygen and bulk defect density of the silicon monocrystal after the heat 
treatment based on the input data. 

2. A method of determining oxygen precipitation behavior in a silicon monocrystal according to Claim 1 , characterized 
in that said computer is programmed based on a model in which oxygen precipitation is expressed as a function of 
the oxygen concentration of the silicon monocrystal and the temperature and time of the heat treatment, and free 
energy produced due to formation of oxide precipitates depends on the impurity concentration of the silicon monoc- 
rystal. 

3. A method of determining oxygen precipitation behavior in a silicon monocrystal according to Claim 2, characterized 
in that said free energy depending on the impurity concentration is strain energy. 

4. A method of determining a process for producing silicon monocrystalline wafers, characterized in that the initial- 
oxygen-concentration dependency, impurity-concentration or resistivity dependency, and crystal -thermal-history 
dependency of an amount of precipitated oxygen and bulk defect density are determined in accordance with the 
method according to any one of Claims 1-3; and an initial oxygen concentration, impurity concentration, resistivity, 
and thermal history of the silicon monocrystal are determined such that wafers will have a desired amount of pre- 
cipitated oxygen and bulk defect density during or after a specific wafer process. 

5. A recording medium on which is stored a program of determining oxygen precipitation behavior in a silicon monoc- 
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rystal through use of a computer, the program instructing the computer to input an initial oxygen concentration of a 
silicon monocrystal, an impurity concentration or resistivity of the silicon monocrystal, conditions of heat treatment 
performed on the silicon monocrystal, and to calculate, based on the input data, an amount of precipitated oxygen 
and bulk defect density of the silicon monocrystal after the heat treatment. 
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